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FOREWORD 


The research reported in this document was motivated by the need 
for a definition of low level wind shear environments for use in 
studies relative to the effect of wind shears on aeronautical 
systems during the landing flight phase. Modeling the wind environ- 
ment has received renewed interest in the last few years because of 
the comnitment on the part of the aeronautical community to develop 
an all-weather automatic landing system. A wide variety of wind 
models has been proposed, each with its own merits and deficiencies, 
depending on the intended application. The model contained herein 
is an excellent contribution to this existing ”stable"of models, and 
was developed for the assessment of the effects of wind shear on 
the landing flight phase of aeronautical systems. The user is 
cautioned against selecting a given model for a design or operational 
problem without examining the other available models. 

This research was conducted by the University of Dayton Research 
Institute for the National Aeronautics and Space Administration, 
George C. Marshall Space Flight Center, Huntsville, Alabama, 
under the technical direction of Dr. George H. Fichtl of the 
Aerospace Environment Division. The support for this research 
was provided by Mr. John Enders of the Aeronautical Operating 
Systems Division, Office of Advanced Research and Technology, 

NASA Headquarters. 
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INTRODUCTION 


An accurate and detailed description of the wind and turbulence 
structure of the planetary boundary layer has many applications. For 
example, a knowledge of wind and turbulence properties can be used to 
determine flow patterns and atmospheric mixing properties for use in 
pollution control. In addition. Structural Engineers must now consider 
the response of large buildings to the wind field; this requires a vertical 
profile of the wind field. The following discussion concerns the aircraft 
landing problem, particularly with respect to STOL aircraft. The two 
components of this problem are (a) the wake effects of buildings around 
STOL ports, and (b) the effects of a sudden change in wind magnitude or 
direction during final approach. A vertical wind profile over and around 
the structure in conjunction with information about wakes behind buildings 
is needed as input in order to define the former effect. A study of the 
latter effect has resulted in the formulation of a structural model of wind 
and turbulence in a statistically stationary and horizontally homogeneous 
boundary layer which has been used as the wind input to simulated air- 
craft landings. By observing the deviation in touchdown point from a 
zero wind condition, the types of profiles that provide hazardous landing 
conditions can be analyzed. 

GENERAL DISCUSSION 

A model has been selected for the mean wind profile and for the 
spectral characterization of turbulence in the statistically stationary and 
horizontally homogeneous surface boundary layer. In the surface boundary 
layer, only (a) the surface conditions, (b) the stability conditions, and (c) 
altitude are considered to influence the wind structure. The horizontal 



shearing stress is constant. The wind direction is considered constant 
with altitude. Inertial, Coriolis, and pressure gradient forces are not 
considered significant in the surface boundary layer. The surface layer, 
in a strict sense, extends to only about ten meters altitude, but the above 
surface layer assumptions are reasonable to an altitude of nearly 100 
meters. 


SUMMARY OF MEAN WIND MODEL 

The mean wind profile that has been selected for the surface 
boundary layer is a function of the surface parameters, u* (surface friction 
velocity), Z Q (roughness length), and d (zero plane displacement); the 
stability parameter Z/L (L is the Monin-Obukov stability length); and the 
height, Z. The model defined by these parameters reasonably approxi- 
mates the true wind profile to an altitude of 100 meters. A definition of 
the parameters u*, Z Q , d, Z/L, and a discussion of rationale pertinent 
to the selection of this model is presehted in a later section. 

The mean wind for the surface boundary layer is described by 
the Logarithmic law 

u = ^~[ in (Z+l2.) + 1(2)] 

where, for neutral stability, 

f (|).0; 

for an unstable condition. 



with the parameter a = 18; 
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and for the stable condition, 

Z Z 

¥( — ) = a — , with the parameter a = 5. 2. 

For the very stable condition, the logarithmic wind law does not 
hold. Under strong inversion conditions, the layers of the atmosphere 
become disconnected and large scale frontal motions are the predominate 
factor in defining the shape of the profile. Figures 5 and 6, shown and dis- 
cussed later, define typical types of profiles that may occur in the very 
stable condition. 


SUMMARY OF TURBULENCE SPECTRA 


Deviations in wind velocity having periods of less than ten minutes 
are considered gusts or turbulence. The statistical properties of turbulence 
are characterized by its spectrum. For the aircraft landing simulation 
problem, the functional representation of the turbulence spectra has been 
chosen as the Dryden spectral function. The parameters of the Dryden 
spectra, tr u , cr v , a w , and L u , L v , L w , have been defined as functions of 
stability, height, and surface conditions. The surface roughness para- 
meter, Z Q , is not sufficient to explain the terrain effect on the components 
of <r and L. Large scale terrain features are important, but insufficient 
data has been collected to categorize their effect. A later section discusses 
the various sources of data which served as a basis for choosing this model. 


The expression for the vertical velocity standard deviation is 


O' 


w 


1. 25u* ( 1 


Z/L 1/4 
S ' 
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where 


and 


S = 1 + 5. 2 J- for L > 0 

S=(l-18~ )“ 1/4 for ~ < 0 
L» Lj — 


The longitudinal and lateral velocity standard deviations, cr^ and cr v , are 
defined in terms of the stability parameter B by Figures 8 and 9, presented 
later, where 


B = 


Ri S 


i n(~)-Y(^-) 
^o u 


The vertical scale length is defined by 

0. 74Z 


L w = 


$ 


whe re 


Z Z 

$ -1+9 — for — > o 

G L L 


and .$ = (1 - 18^-)" 1/4 --L for L < 0 

£ L Li Li 


The longitudinal and lateral scale lengths are derived from the local 
isotropic turbulence relationship 


L ZL ZL 

U W V 



The and defined by this relationship show generally good agree- 
ment with independently measured experimental data. 


SURFACE BOUNDARY LAYER 

The surface boundary layer is defined that region of the atmos- 
phere where the shearing stress is constant. The region extends, in a 
strict sense, from the surface to only about ten meters altitude, but is 
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considered to extend to nearly 100 meters. In the ten-meter to 100-meter 
layer, the decrease in stress with altitude is so small as to make the 
assumption of constant shear stress acceptable for engineering applica- 
tions. In the surface boundary layer only the surface conditions, the 
stability condition, and altitude affect the wind and turbulence structure. 

The following sections discuss the parameters used to describe the sur- 
face and stability conditions. 

Surface Conditions 

The surface parameters used in defining both the mean wind 
velocity profile and the spectral properties of turbulence are the surface 
roughness length (Z 0 ), the zero wind reference level (d), and the sur- 
face friction velocity (u*). The roughness length Z Q is a parameter used 
to characterize the gross features of the terrain. It is most reliably 
estimated for uniform type terrains. When the terrain is nonuniform, 
the dissimilar features influence the wind profile in a manner that depends 
upon the distance from these features. Thus, a gross roughness para- 
meter is not appropriate for estimating a wind profile if the profile 
reflects the effects of nonuniform terrain elements. Fichtl [ 1 ] gives 
some typical values of roughness length for various types of uniform 
terrains. (See Table 1. ) 

The zero plane reference level, d, is the altitude at which, for 
a given roughness length, the wind velocity extrapolates to zero. (See 
Figure 1. ) For example, over forest type terrain where the average 
height of the trees may be 15 meters, the roughness length is Z Q - 
0. 5 meter. The profile above the trees can be represented as a logarithmic 
function of altitude. The zero plane displacement is determined by 
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TABLE I* 


TYPICAL VALUES OF SURFACE ROUGHNESS LENGTH (Z Q ) 
FOR VARIOUS TYPES OF SURFACES 


Type of Surface 

N 

o 

B 

Mud flats, ice 
Smooth sea 
Sand 

Snow surface 

Mown grass (~ 0.01 m) 

Low grass, steppe 

Fallow field 

High grass 

Palmetto 

Suburbia 

City 

io " 5 - 3-10' 5 
S-IO -4 - 3- lo ” 4 
lo -4 - io " 3 
10 " 3 - e-icf 3 
io -3 - io ' 2 

10 ~ 2 - 4’ 10 2 

-2 -2 
2*10 - 3*10 

4-10~ 2 - 10 _1 

10' 1 - 3-iO -1 

1 - 2 

1 - 4 


*This table was taken from George H. Fichtl, Reference 1. 
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Figure i. Logarithmic Wind Profile Above a Forest. 


7 






extrapolating the wind profile above the tree tops to the altitude at which 
the wind equals zero. In moderate wind fields, the zero plane displace- 
ment may be the approximate height of the roughness elements, in the 
above example, approximately 15 meters. In strong wind fields, the zero 
plane displacement tends to decrease while in light wind fields, it tends 
to increase. Thus, for an average guess, d can be taken as the average 
height of the roughness elements. 


r — ! 1/2 

The surface frictional velocity u* is defined as u* = [ -uw J , 
where w and u denote the vertical and longitudinal components of the 
fluctuating wind velocity, and overbar denotes a time average over approxi- 
mately a ten-minute period. The fluctuating velocity is distinguished from 
the mean velocity by its duration. In practice, the mean wind is generally 
averaged over periods approximately equal to ten minutes and the fluctu- 
ating velocity is computed as deviations from the mean. Since a gap in 
spectral energy exists between approximately ten minutes and several 
hours, there should be little difference between an average computed over 
ten minutes and one over a one-hour period. [2] Thus, it is reasonable 
to define a fluctuating component of velocity as having a period less than 
ten minutes. 


The frictional velocity is a scaling parameter for the mean wind 

profile. That is, the mean wind profile increases in direct proportion to 

the surface friction velocity. This can be seen as follows: The corre - 

lation coefficient, p = uw /cr <r , is assumed to be constant throughout 

uw u w 

the surface boundary layer (typically, p = -0. 3). Thus, uw = Dcr a* , 

uw u w 

D= -0 o 3, a constant for this discussion. In neutral stability cr and cr 

u w 

are known to vary proportional to the mean wind velocity. That is, 

cr - Au and cr = Bu . Thus , u* = [ -uw ] ^ ^ = Cu, whe re C = \ DAbT ^ . 
u w L J 9 J 


8 



Thus, in neutral stability, the mean wind is directly proportional to the 
surface friction velocity. In nonneutral stability, similarity theory again 
requires that u be directly proportional to u*. 


Stability Parameters 


The mean wind and turbulence in the surface layer depend upon the 
temperature profile. The temperature lapse rate, i. e. , the decrease in 
temperature with altitude, measures the effect of the temperature profile 
on the wind velocity structure. If the decrease of temperature with altitude 
is greater than that associated with a homoentropic atmosphere (adiabatic 
lapse rate) then the temperature profile adds kinetic energy to the atmos- 
phere via positive buoyancy forces. If, on the other hand, the temperature 
lapse rate is less than the adiabatic lapse rate, a downward restoring force 
absorbs energy from the atmosphere. In this case, the atmosphere is said 
to be stable. If the temperature lapse rate is approximately adiabatic, the 
atmosphere is neutral. Under strong inversion conditions, the damping 
effect becomes so prominent as to convert turbulent eddies into a laminar 
flow. This division between laminar and turbulent flow is generally defined 
by the critical Richardson number, Ri . For Ri > Ri , turbulence can no 
longer exist. This critical value is not accurately defined and may even 
be dependent on surface condition. An approximate value for the critical 
gradient Richardson number is Ri=0. 20 [2]. The expression for gradient 
Richardson number in terms of lapse rate and the adiabatic lapse rate is 


Ri 


JL 

T 


(£ - V) 

<— ) 2 

'*Z ' 


(1) 
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where 


g = gravitational acceleration, 

F = adiabatic lapse rate = g/Cp , 

■y = temperature lapse rate = -^T/3Z, 

T = temperature, 

Cp = specific heat of dry air at constant pressure, 
u = mean wind velocity, and 
Z = altitude. 

Equation (1) shows that for neutral stability Ri~0; for an unstable atmos- 
phere Ri < 0; and for a stable atmosphere Ri > 0. The stable atmosphere 
is further subdivided by the critical Richardson number into stable for 
Ri < Ri^ and very stable for Ri>Ri^. In the very stable condition, when 
only laminar flow exists, the atmosphere becomes disconnected and the 
shape of the wind profile no longer reflects only the surface conditions. 
Large scale effects such as a warm front overrunning a cold front become 
the dominant features in shaping the profile. Physically, the Richardson 
number represents the ratio of the thermal to mechanical production of 
turbulence. For a statically stable atmosphere, T >■ y, thermal damping 
occurs, and the Richardson number is positive. For an unstable atmos- 
phere, y >F, thermal convection adds energy to the atmosphere, and the 
Richardson number is negative. For a neutral atmosphere, there is 
essentially no thermal energy transferred and the Richardson number is 
approximately zero. The Richardson number has become an important 
and universally accepted measure of atmospheric stability. In defining 
wind profiles, however, the Richardson number is not easily manipulated 
because of its variability with altitude. A different parameter, Z/L, 
where L is the Monin-Obukov stability length, is a more convenient 
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stability parameter for characterizing wind profiles. The Monin-Obukov 
stability length can be considered independent of altitude in the surface 
layer and is defined by 

T __ u*3Cpp0 
L= kgH 

where 0 = potential temperature, 

H = vertical heat flux, 

P = air density, and 
k = von Karman constant. 


The parameter Z/L can be related to the Richardson number through 
dimensional analysis arguments and experiments show the relationship 
to be given by 


Z/L = Ri for Ri< 0 (Businger's hypothesis) 

and 

Z/L = Ri/1-5. 2 Ri for 0 < Ri < Ri . 

c 


The second expression above implies Z/L -+co as Ri -* 1/5.2. 

This conforms to a critical Richardson number of approximately 0. 20. 


MEAN WIND PROFILE 


In the surface boundary layer the logarithmic wind law as derived 
from the Monin-Obukov similarity theory and experiments for neutral 
(Ri=0), stable (0< Ri< Ri c ), and unstable (Ri < 0) conditions is given by 


u 


u* 

k" 


U n(- 


Z'-d+Z, 


) +M 


Z' - d 
L 


)] , 


(2) 
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where the surface parameters u*, Z Q , and d, and the stability parameter, 
Z/L have been previously discussed. Z 1 is the height above the surface. 
The von Karman constant, k, has numerical value of approximately +0.4. 
Equation (2) can be written in the simplified form 


u 


u* 

IT 


Z + z 

[1 n (-—> + *<!>] 


(3) 


where Z is the altitude above the zero plane reference level. Equation (2) 
or (3) is useful only for fully turbulent boundary layers and does not pro- 
vide a valid representation of the wind profile under very stable conditions. 


The function V (Z/L) is an empirically derived universal function 
of Z/L. For neutral stability when Z/L = 0, v (0) = 0, so that the log 
profile is valid, i„ e. , 


u 



Z + Z 

o 

~Z 

o 


) 


Figure 2 shows neutral wind profiles for various values of u* and 
Z . The expression for v (Z/L) in unstable and stable air has been derived 
by many authors (Webb, Panofsky et al. , McVehil, Fichtl and McVehil, 
Wyngaard and Cote, Businger et al. , and others ). [2, 3, 4, 5, 6, and 7] 

For the unstable conditions ¥ (Z/L) has been calculated by fitting 

k Z \1 

the nondimensional shear S = — to an equation of the form 

u* oZ 

S = (1 - a Ri)‘ 1/4 (4) 


where Ri = Z/L and a is a fitting parameter. 
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Altitude 
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0 10 20 30 _ 40 50 60 70 ft/sec 
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Figure 2. Neutral Wind Profiles. 
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Peterson and Panofsky [ 3 ] suggest a value of a = 18; Businger 
et al. suggest a = 15 withavon Karman k of 0. 35. Integration of the non- 
dimensional shear (Equation 4) defines the mean wind as 


u 


J, Z t Z rj 


where 




Figure 3 shows unstable wind profiles for various values of Z q . u*. and 
L. 


For a stable atmosphere, the nondimensional shear has been des- 
cribed by the expression 


S = 


kZ 9u , Z 
u* BZ “ + a L 


(5) 


Integration of Equation (5) with respect to altitude defines the mean wind 
profile for stable air as 


u 


Z + Z 

= TT 0 “ z ~ 




O 

Peterson and Panofsky suggest a value for a of 10; Businger et al. 
suggest a = 4. 5 to a = 5. 0 with k = 0. 35; McVehil suggests a= 7. 0. Webb 
suggests a = 5. 2 with individual observation showing a standard deviation 
of 30 percent from the mean value of 5. 2. There is, in essence, some 
doubt as to whether the nondimensional shear is dependent only upon Z/L. 
For our wind shear model, the value of Webb, a = 5. 2, has been selected. 
Figure 4 shows plots of stable wind profiles for various values of u*, Z , 
and L. 
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For the very stable conditions no analytic expression has been 
found to represent the wind profile. The shape of the profile is largely 
determined by large scale motions of a new air current overrunning a 
static front. Since the flow is laminar, very little mixing occurs and 
the overriding air remains separated from the underlying air. This 
condition also leads to the possibility of a large change in wind direction 
in the surface layer. Changes of wind direction in excess of 45 degrees 
are not uncommon in very stable air. [8] A few examples of the types 
of wind profiles that may occur are given in Figures 5 and 6. Figure 5 
shows a calm below the interface level Z^. Above the interface the wind 
is considered constant of magnitude u, depending upon the velocity of 
the overriding air. u can be assumed constant above because the 
boundary effect produced by the earth's surface is not relevant and the 
i boundary effect between nonmixing air masses of different velocities is 
negligible. Figure 5 also allows for a variation in the altitude Z where 
the interface occurs. Figure 6 shows a light wind condition which obeys 
the logarithmic wind law below the interface with a constant wind above 
the interface . This situation may occur when the atmosphere is in a 
neutral or stable condition below the interface. Large changes in wind 
direction are likely to occur with this type profile. 

TURBULENCE SPECTRA 

Measurements of deviations from the ten-minute to one -hour 
averaged mean wind are classified as turbulence or gusts. These gust 
measurements cannot be characterized in a deterministic sense and 
hence must be described statistically. The turbulence structure of the 
atmosphere is generally specified by its power spectra. Studies of 
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Figure 5. Very Stable Wind Profiles. 



Figure 6. Very Stable Wind Profiles. 
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measured spectra have been fit to functions so that a functional repre- 
sentation of the spectra is possible. Some of the more prominent spectral 
density functions are the Dryden, von Karman, Busch-Panofsky, and 
Fichtl-McVehil spectral functions. The Dryden and von Karman functions 
are defined in terms of the parameters 0 - and L. Sigma is the standard 
deviation of the gust velocity and L. is the integral scale length of tur- 
bulence defined as the integral of the correlation function. The Busch- 
Panofsky function is defined by the parameters <r and f , where f 
is the nondimensional frequency at which the spectrum takes on its 

maximum value, f is a function of stability only. The Fichtl- 

max 

McVehil spectrum is a modification of the Busch-Panofsky spectrum 
which contains an additional parameter which is dependent upon height. 
Much of the data that has been collected has been fit by the Dryden and 
von Karman functions. From this data variation in cr and L. with terrain, 
stability, and height has been analyzed. On only a relatively small per- 
centage of the data has the variation in f with stability been examined. 

max 

The von Karman spectrum differs from the Dryden spectrum primarily 
in the inertial subrange region. In this region local isotropy holds and 
there is essentially no production or dissipation of energy. The spectrum 
is proportional to the wave number, K, to the -5/3 power. Experimental 
data has largely substantiated the local isotropy and -5/3 slope hypothesis. 

The von Karman function has the desirable property of being pro- 
-5/3 

portional to K in the inertial subrange. The Dryden function is pro- 
-2 

portional to K in the inertial subrange. Thus, from this standpoint, 
the von Karman is the more desirable spectral density function. 

For the aircraft landing problem, for which use of this model is 
intended, a digital simulation of turbulence is required. The von Karman 
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function, being an irrational function, is more difficult to simulate than 
the rational Dryden function. Furthermore, it has been shown by Flinn 
[9] that in observing aircraft response little difference can be seen 
between a Dryden and a von Karman input spectrum. Consequently, for 
the landing simulation problem, the Dryden spectrum has been chosen. 

The parameters of the Dryden spectrum cr and L are functions of 

stability, surface conditions, and altitude. The standard deviation of a 

gust velocity time history from a running ten-minute to a one -hour mean 

is a common method used to measure cr. Consequently, the determination 

of cr is independent of the spectral density function being used (e. g. , von 

Karman, Dryden, etc. ). The same is not always true for the scale length 

parameter, L. The scale length L is often calculated by observing the 

wave number, K , at which the spectrum is a maximum and then re- 
max 

lating K to L through the relationship 
max 

L = 0. 146/K 
u max 

and (6) 

L =0. 106/K = L for the von Karman spectrum, 

v max w 

and 

L = 0. 159/K 
u max 

and (7) 

L = 0. 1 1 7/K = L for the Dryden spectrum. [ 1 0 ] 

v max w 

Since Equations (6) and (7) are not identical, it follows that when 
fitting an experimental spectrum with the von Karman and Dryden functions, 
slightly different values of L would result. Percentagewise, however, 
the difference between the L values is small compared to the uncertainty 
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in defining K^. Thus, for practical considerations, experimental 
values of L. can be analyzed without regard to which spectral function was 
used to obtain L. If L values were obtained by an entirely different pro- 
cedure such as by direct integration of the correlation function,* they would 
be independent of the type of spectral function. 

The following sections examine the stability, terrain, and height 
dependence of cr’s and L's obtained from experimental measurements. As 
concluded in the above paragraph, the c r and L relationships are valid, for 
use with both the Dryden and von Karman models. 

Standard Deviation of Gust Velocities 

Vertical Standard Deviation - Similarity theory predicts that the 
standard deviation of the vertical component of gust velocity normalized 
by the friction velocity depends only upon the stability parameter Z/L, 
i. e. , 

~ = f (Z/L). 

U* 

The function f as derived by Monin is 

£ = ct i -^>] i/4 (si 

where S is the nondimensional shear defined by Equations (4) and (5). 
Prasad and Panofsky [3] have shown Equation (8) to provide a good fit 
to experimental data when c = 1. 25. This value of c implies that for the 
neutral stability condition 


21 



= 1. 25, 


c r 

w 
u* 

a value that is in good agreement with what several other experimenters 
have observed in neutral stability. [3, 11, and 12] Other expressions 
can be found to represent the variation of c with stability, most of which 
are in general agreement with Equation (8). Thus, Equation (8) will be 
used to define the standard deviation of the vertical gust velocity component. 
Figure 7 shows a plot of Equation (8). 

Longitudinal Standard Deviation - The longitudinal component of 

gust, in contrast to the vertical component, has not been found to obey 

similarity theory. In particular, the ratio cr /u* is not a function of 

stability only, but is also dependent upon large scale (hills, trees, 

buildings, etc. ) surface terrain features. The surface roughness length, 

Z , does not appear sufficient to explain the variation in the ratio cr /u* 
o u 

from one place to another in a given stability condition. The author is 
not aware of any study that has explained the terrain effect on cr^/u*. In 
neutral stability various measurements of the ratio have provided values 
between 2. 1 and 2. 9. [ 13, 14] The value of cr /u* = 2. 5 will be taken 

as an average value for this model in neutral stability. 

The variation of the ratio <r / u* with stability and height has been 
analyzed by Prasad and Panofsky [3] by introducing the stability para- 
meter B defined as 

2 

B = (fH-gzHr-v) • 

B is related to Z/L, Z^, and Z through the expression 
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(9) 


B = 


RiS 


i n (|- o )-^) 


where the parameters are as defined previously. Prasad and Panofsky 

show the relationship between c r /u* and B in Figure 8. In neutral sta- 

u 

bility, Z/L = 0, B = 0, and cr /u* = 2. 5. The variation of cr /u* with B 

u u 

is small, especially when compared with the lateral component. Although 

not immediately obvious by examination of Figure 8, other parameters 

constant, the variation of cr /u* with Z/L. is also small. Other investi- 

u 

gations have found the same type of variation in %r fu* with stability and 
height as Prasad and Panofsky. Since Prasad and Panofsky's work en- 
compasses a wide variation in stability conditions and provides mathe- 
matical expressions for calculating o r , it has been chosen for inclusion 
in this model. 

Lateral Standard Deviation - The lateral component standard 

deviation of gust velocity, like the longitudinal component, does not obey 

similarity theory. The ratio of o* /u* in neutral stability has been found 

v 

to vary from 1. 3 to 2. 6, depending not only upon surface roughness, Z^, 

but apparently also upon the large scale terrain features. [ 15] The 

author knows of no mathematical expression that relates cr /u* to Z 

L v o 

and large scale roughness in neutral stability. Hence, an approximate 

value of a* /u* = 2. 0 in neutral stability will be used. The variation of 
v 

cr / u* with stability is much greater than for the longitudinal components. 
Prasad and Panofsky show the variation of o* / u* with B in Figure 9, 

At Z/L = 0 , B - 0, and <r /u* = 2. 0. Percentagewise, the variation of 
cr /u* with stability is large, particularly in the unstable condition. 

Other investigators agree in general with Prasad and Panofsky 1 s modeL[5] 
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Figure 8. The Ratio cr / u* as a Function of B at Many Sites. 
(From Pano¥sky et al. [3]) 
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Figure 9. The ratio cr /u* as a Function of B at Many Sites. 
(From Panofsky et al. [3]) 
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Turbulence Spectra Scale Length 


Vertical, longitudinal, and lateral components of scale length 
can be measured independently from experimental data. Physical laws, 
however, may demand an interrelationship between the various com- 
ponents. For example, in the free atmosphere, above the boundary layer, 

it is reasonable to apply the isotropic constraints cr = cr = cr and L = 

r u v w u 

2L^ = 2L^. In the surface boundary layer, total isotropy is certainly 
not valid, but local isotropy, or isotropy at high frequencies, is a tenable 
assumption. Total isotropy implies that the relationship 


$ 

w 


(K) 


is K " x u 

2 u ” 2 3K 


= $ (K) 

v 


( 10 ) 


holds for all wave numbers K. For local isotropy, we will only require 
that Equation (10) hold for large wave numbers; that is, as K +00 • The 
longitudinal, vertical, and lateral spectral functions for the Dryden 
spectrum are 

4a- L 

$ = , UD 


u 1 + (2 ttL K)4 
u 


and 


$ 

w 


I 

V 


4c r ^L [ 1+3(4ttL K) 2 ] 
w w w 


[ 1 + (4tt L w K) 2 ] 2 
4cr Z L [ 1 + 3( 4ttL K) Z ] 

V V V 

[ 1 +(4tt L K) 2 ] 2 

V 


( 12 ) 

(13) 


Substituting Equations (11), (12), and (13) into Equation (10) and letting 
K *oo yields the relationship for local isotropy 


L 

u 


2 


cr 


u 


2L 

w 


2 

cr 

w 


2L 

v 

2 


or 

v 


(14) 
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Having already described cr^, and cr for the surface boundary layer, 

it is only necessary to specify one component of L and then apply 
Equation (14) to define the other two components of scale length. 

Vertical Scale Length - The vertical component of scale length 
can best be described independent of the local isotropy constraint for two 
reasons: (a) a large amount of experimental data exists concerning the 

vertical scale length, and (b) the conformance of the vertical component 
of turbulence to similarity theory implies a linear relationship of vertical 
scale length to height. 

The vertical scale component is affected by the large scale terrain 
features. The surface roughness length may have some influence on L^, 
but the nonuniformities and large scale features predominate. For neutral 
air over relatively flat terrain, a number of experimental results are 
available. Figure 10, extracted from Teunissen [ 10 ] , shows results 
from several investigators. Teunissen summarizes their results and 
finds the relationship 

L = 0.4 Z (15) 

w 

to provide a reasonable fit to the combined set of data. Sufficient results 
are not available to determine how the proportionality constant of Equation 
(15) varies over different types of terrain. 

The variation of the vertical scale length with stability has been 
examined by Busch and Panofsky. [3] Busch and Panofsky find that the 
nondimensional frequency at which the normalized spectrum takes on its 
maximum value is related to stability by the expression 

f = K Z = 0.32$ (Z/L) (16) 

max max e 
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where 


(Z/L) is the dimensionless dissipation rate. 

For the Dryden spectrum, the vertical scale length, L^, is related 
to the wave number where the normalized spectrum is a maximum by 
the equation 

L = 0. 117/K 
w max 

or T 0. 117/Z ‘(17) 

w f 

max 

A substitution of Equation (16) into Equaton (17) gives 

T 0. 37Z 

W " $E 

An analysis of several sets of data by Busch and Panofsky has provided 
the solid curve shown in Figure 11 (obtained from Reference 3) as repre- 
senting the nondimens ional dissipation rate ? . Two hypotheses have 

been suggested that postulate a balancing of the dissipation rate with 
other forces. They are 

$ = S (18a) 

S 

and _ 

$ = S - £ ; (18b) 

e L 


that is, 

(a) dissipation rate equals mechanical energy production 

and 

(b) dissipation rate equals mechanical plus buoyant production. 
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ord and Cramer (1966), corrected 




Figure 11 also compares Busch and Panofsky's curve to Equation 

(18a) and (18b) for S defined by Equations (4) and (5). On the stable side, 

Z 

neither assumption § = S or § = S - — fits the data. Busch and Panofksy's 

solid line is derived from McVehil f s[4] nondimensional shear expression 
Z 

of S = 1 + 10— and Equation (18b). This expression for the nondimensional 

jLj 

shear is inconsistent with recent findings by Webb[Z] and Busingeret al. 
[6] Consequently, on the stable side we will choose Busch and Panofsky's 
function 

$ e = 1 + 9 | (19) 


for the dissipation rate because it empirically fits the data, while rejecting 
both balancing assumptions (18a) and (18b). 


On the unstable side the dissipation rate appears balanced by the 
sum of the mechanical and bouyant energy production, at least to Z/L = 
-0.4. In strong instability sufficient data is not available for analysis. 
The function 


$ 


e 



(1-185)- 1 ' 4 . 


Z 

L 


( 20 ) 


is used in the model to expressthe dissipation rate for the unstable condition. 


The expression for the vertical scale length is obtained by sub- 
stituting Equations (19) and (20) into Equation (17). This gives 


and 


L 

w 


0. 37 

1+9 i 



( 21 ) 


L. 

w 


0. 37 

1174 " 


d-1 8 !> 


Z 

L 
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In neutral stability, when Z/L = 0, Equations (21) reduce to = 
0. 37Z which is in good agreement with the proportionality constant of 0. 4 
recommended by Teunissen. 

Longitudinal Scale Length - The longitudinal scale length is de- 
rived from the vertical and longitudinal standard deviations and the 
vertical scale length using the local isotropy relationship 


L 

u 


2 


O' 

u 


2L 

w 



(14) 


A comparison can be made in neutral stability to determine how well the 
values of L^ derived from Equation (14) compare with experimental mea- 
surements. In neutral stability the solution to Equation (14) for L gives 

L = 2. 96Z 
u 

Figure 12, from Teunissen [ 10 ] , compares the line L = 2. 96Z with 
various experimental results. Although the line could not be considered 
an excellent fit, it does reasonably express the variation of L^ with 
stability, particulary above three-meter altitude. 


Lateral Scale Length - As in the case of the longitudinal scale 
length, the lateral scale length is derived from the local isotropic rela- 
tionship 

L 2L 


v w 
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In neutral stability, the solution to the above equation is 


L v = 1. 89Z. (22) 

Figure 13 compares Equation (22) to the relationship found by Fichtl and 
McVehil in neutral air between 18 and 150 meters. The slope of the two 
curves disagrees considerably even though they pass through the same 
value at approximately 40 meters. Little other data is available from 
which to draw conclusions. 
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